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Abstract 

We present results from extensive numerical integration of the KPZ equation in 1 + 1 dimensions 
aimed to check the long-time behavior of the dynamical structure factor of that system. Over a 
number of decades in the size of the structure factor we confirm scaling and stretched exponential 
decay. We also give an analytic expression that yields a very good approximation to the numerical 
data. Our result clearly favors stretched exponential decay over recent results claiming to yield the 
exact time dependent structure factor of the 1 + 1 dimensional KPZ system. We suggest a possible 
solution to that contradiction. 
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Many interesting dynamical phenomena in condensed matter physics are described in 
terms of non-hnear field equations driven by noise. A long list of examples includes turbu- 
lence, critical dynamics, the dynamics of interacting polymers, ballistic deposition (as well 



as other growth models) etc. The Kardar- 



growing surface under ballistic deposition 



m 



arisi-Zhang (KPZ) equation that describes a 



is such a model. This equation formulated 



in terms of a height function h (r , t) driven by external noise is given by 

dh(r,t) „9, A,„,,2 N /,N 

' = i^y^h + - {Vhf + ?7 (f, t) , (1) 

where v is a. diffusion constant, A is the coupling constant (that controls the sticking rate 
of the deposited material), and t] (r, t) is a noise term driving the equation that models the 
randomness of the falling material. The noise is usually chosen to be Gaussian, with zero 
mean and second moment 

{r, (f, t) 77 (r^, t')) = 20,5" (f - r') 5 [t - t') , (2) 

where d is the substrate dimension and -Do specifies the amplitude. 

The KPZ equation has been suggested 17 years ago I] as an extension of the linear 
Edawards- Wilkinson equation j^, so that a lot of research has been done on the statistical 
properties of the surfaces that this equation grows. An extensive review of this work can 
be found in refs. It is well known that KPZ surfaces are self-affine, and are 

well described by two scaling exponents, namely the roughness exponent a and the dynamic 
exponent z. It turns out that in the KPZ system these two exponents are not independent. 
Due to symmetry of eq. with respect to infinitesimal tilting 1^ the famous scaling relation 
a + z = 2 is established. Furthermore, for the special case when d = 1, the existence of a 
fiuctuation-dissipation theorem gives the exact result a = 1/2 and z = 3/2. 

In recent years there has been a growing interest in the dynamical properties of the 
KPZ system. A question of great interest regards the long-time behavior of the dynamical 
structure factor $q (t) = {hg (0) (t))g, where hg (t) is the Fourier transform of the height 
function h{f,t), and (■ ■ ■)g denotes steady-state averaging over the noise. Notice that by 
definition $g (0) = (/)g, where (pg is the static two-point function. 

Using a self-consistent approach, Schwartz and Edwards were able to predict a stretched 
exponential decay for ^g{t) 0, 0]. Regarding the KPZ system, they found the following 
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long-time asymptotic behavior 



d-l 



%{t)r^c<P,{jqt'/^) ' e-^^*"\ (3) 

where c and 7 are dimensionless constants (not necessarily 1), d is the dimension, and z 
is the dynamic exponent. The same asymptotic behavior was also predicted analytically 
by Colaiori and Moore using a mode-coupling approach. Later, Colaiori and Moore 
solved numerically the mode-coupling equations in one dimension and confirmed the 
asymptotic analysis for the long-time behavior. Surprisingly, they also found that $y (t) 
decays to zero in an oscillatory manner - a fact that was not revealed by the analytical tools. 

It should be stressed however that the above describe only approximations to a solution 
of the KPZ equation. Even an exact solution of either the Self-Consistent Expansion (SCE) 
equation or the mode-coupling equation would provide only an approximation for the real, 
time-dependant structure factor, $q (t), of the KPZ equation. Indeed, a more recent publi- 
cation by Prahofer and Spohn suggests that the envelope of the one-dimensional structure 
factor decays exponentially rather than as a stretched exponential . They claim an exact 
solution for the time- dependant structure factor of another model in the universality class of 
the KPZ system, namely the polynuclear growth model. The actual solution for the struc- 
ture factor follows a number of well defined steps involving some direct though complicated 
numerical calculations. Those calculations are reported to be performed with extremely 
high precision that seems to ensure that the final solution for the structure factor is not 
affected by inaccuracies in the numerical procedure. Although the solution they present is 
numerical it is rather obvious that the envelope decay is exponential rather than stretched 
exponential. 

This discrepancy motivated us to check the above results directly on the KPZ equation. 
In this work we find numerical support for the existence of stretched exponential relaxations 
in the KPZ system in 1 + 1 dimensions in contradiction to ref. 3|- We also get direct 
evidence for the oscillatory behavior of As a by product, we were able to verify 

the predicted short-time behavior of the dynamical structure factor (given in ref. for 
example), and the validity of the scaling hypothesis for small g's. 

We discretized the KPZ equation on a one dimensional lattice, with lattice constant 
Ax, and time difference At, 

h{x,t + At) = h{x,t) + 



H 2 y2{^[hi^ + t) - 2h {x, t) + h{x- Ax, t)] 

(Ax) 1^ 

+ ^[h{x + Ax, t)-h{x- Ax, + a (UAt)^^^ t] (t) , (4) 

where = 2Dq/ Ax and the random numbers t] {t) are uniformly distributed between —1/2 
and 1/2. In this work we used L = 1024, At = 0.05, Ax = 1, = 1, A = 4 and Dq = 1/6. 
After reaching steady state, at each time step we Fourier transformed the discrete height 
function, and obtained hg (t) for q = qQ,2qQ, . . ., where go = 27r/1024. Then we calculated 
the two-point function {hg (0) h^g (t))^ by averaging over all pairs with a time difference t, 
for q = 30qo, 60qo, 120go and 240go- 

First, the scaling hypothesis was numerically verified, i.e. In Fig. QJwe plot f (ugt) = 
$g {t)/(f)g as a function of Ugt = Bq^t with z = 3/2 for q = 30go, 1 = 60go and q = 120go- 
The plot indicates good scaling. 




FIG. 1: A log plot of the scaling function / {tOgt) for various small g's. 

However, we found that the scaling hypothesis breaks down when taking q = 240go (see 
Fig. 121) which means that the scaling, which is supposed to be correct for small g's, is correct 
only up to g's that are of order of ~ 10% of the largest g in the system. 

Fig. HI indicates good scaling at least up to 120go. Therefore, we invested most of the 
computational effort in this Fourier component, since eq. Q indicates that the larger g 
we take, the faster computational time evolution we get. For g = 120go we took 5 ■ 10^° 
integration time steps of At. Taking this component we found clear evidence for oscillatory 
behavior as shown in Fig. IHl 

The error estimate in Fig. El was obtained using the imaginary part of {hg (0) h_g {t))g 
after averaging. Note that {hg (0) h^g (t)) ^ should be real, due to averaging, while each 
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FIG. 2: A log plot of the scaling function / (ujqt) for q = 120go and q = 120go- 
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FIG. 3: A log plot of / (ujqt) for q = 120go with an error estimate. 



contribution of the form hq (0) h^q {t) is certainly not real. This means that the imaginary 
part, which should vanish eventually, is a sensible error estimate. An independent argument 
for the estimation of the error, which yields the same order of magnitude runs as follows: 
consider the case where there is no correlation at all and estimate the "apparent correlation" 
due to the finite sample. The total number of time steps is = 5 ■ 10^°. Therefore, this is 
also the number of pairs {hg, h_q) separated by a time t = nAt with n <^ N. The measured 
correlation is 

(5) 



{hq (0) h.q (t)) = J^Y.K (^A^) h.q [(z + u) At]. 

In the absence of correlation, the sum on the right hand side of (0) is a sum of N random 
variables. Therefore, the size of the apparent correlation is of order A^~^/^, which is of the 
order of e~^^. 

Using this result we extracted the short-time behavior of the scaling function. In ref. 
Colaiori and Moore predict (t) oc 0^ 1 — {uOq tf/']. Now, in one dimension F = 2 (F is 
related to the roughness exponent a via V = d + 2a) and 2; = 3/2, so that for a specific q 
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we get / (ujqt) oc 1 — (ujqt) ' . This prediction was indeed verified by our data as shown in 
Fig. H 




FIG. 4: A log plot of / (iOqt) for small iOqfs. 



The error estimates suggest that the data may be useful in the range < Uqt < 55. Since 
we wanted to check a stretched exponential decay, we multiplied our numerical fiojqt) by 
exp [looqtf/^. We chose 7 = 0.93 so as to render the resulting function to be oscillating in 
that region. Motivated by the predictions of ref. Q| and our previous verification of it we 
present the resulting function as a function of {uOqtY/'^ in Fig. 
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FIG. 5: A simple fit is |cos (3 (iOqt) 



n4/3- 



+ D sin 



/3 (wgt)"^/^] } with P = 0.034 and D = 32.7. 



4/3- 



D sin 



(3 [ujqtf^] } for the 



We were thus led to try the fit e'^K*)"'" {cos [(3 {uqt) 
full scaling function f{uqt). This has the right small ujqt behavior as well as the envelope 
stretched exponential decay. 

The fit is very good, especially when keeping in mind that this fit involves more than 6 
orders of magnitude in the size of the scaling function. An attempt to replace the 2/3 in the 
stretched exponential (in the ansatz) with an arbitrary exponent that will be determined 
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FIG. 6: A fit of / (ujgt) using e-^^^^*)"^' {cos [/? (tOgt) 



si/3 



+ D sin 



p (wgt)''/^] } witli 7 = 0.93, 



P = 0.034 and D = 32.7. Botli curves practically join over most of the time range. 

by the fitting procedure yields a very close value (~ 0.65). Actually, we also tried to fit an 
exponential decay (rather than a stretched exponential one), according to the finding of ref. 



12| . that gave a fit that was definitely worse (Fig. [7j). 




FIG. 7: A fit of / (wgt) using e'^^"^'*) {cos [/? (cogt)^^^] + D sin [/? {tOqi f^^^ }. The best fit of that 
form is obtained for 7 = 0.21, /? = 0.326 and D = 2.6. 

At this point we are facing an interesting contradiction. First, we have results following 
from four independent approaches: (1) Analytical asymptotic study of the self-consistent 
approximation. (2) Analytical asymptotic study of the mode-coupling approximation. (3) 
Numerical solution of the mode-coupling equations. (4) The present direct numerical inte- 
gration of the KPZ equation. All these independent approaches yield the same stretched 
exponential decay. On the other hand, we have a publication jjl^ claiming to be exact that 
yields a result that is quite different from the results obtained from all the above methods. 



It is possible that there is a flaw in the study presented in ref. 



121. If there is, we have 



certainly not found any. We would like to suggest here another solution to the problem. Ref. 



12| does not consider directly the KPZ equation but rather the polynuclear growth model 
that was shown to be equivalent to the directed polymer problem with specific boundary 
conditions, and thus in the same universality class of the KPZ system. However, the point 
is that two models that are in the same universality class must have the same exponents 
but not necessarily the same scaling functions. In fact, it is possible to construct families of 
exactly solvable models where all the members of the family are characterized by the same 
exponents yet have radically different time dependant structure factors. We will not go into 
this here but the interested reader could find the relevant ideas, although presented in a 
different context, in ref. Q|. Our suggestion for solving the puzzle is thus that Prahofer 
and Spohn obtain the correct decay for the polynuclear growth model, which is in itself 
a most impressive feat, but this is not the decay of the KPZ structure factor. 

To summarize, using extensive numerical integration of the KPZ equation in 1 + 1 di- 
mensions, this work gives clear support for the scaling hypothesis (i.e. the fact that the 
scaling function / {oOqt) is the same for any q, at least for small g's), verifies the short-time 
behavior given in ref. (i.e. $g (t) oc 0g 1 — {uqt)^^^ ) and establishes the oscillatory 
decay of the dynamical structure factor to zero (as suggested in [ill). In addition, we show 
that the stretched exponential describes the decay of the structure factor over six orders of 
magnitude in its size. This implies that the KPZ problem is likely to be a respected member 
of the family of systems that exhibit slow relaxations - thus opening the door for mutual 
infiuence between the community of surface growth and that of slow dynamics. 

The present work consumed quite a few months of CPU time on a Pentium 4 machine. 
In spite of that, we believe that the results presented here should motivate a much heavier 
numerical effort to deal with a region with larger Ugfs for the one-dimensional system and 
hopefully with the two dimensional system. 
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